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ABSTRACT 

The microwave sky shows unexpected features at the largest angular scales, among 
them the alignments of the dipole, quadrupole and octopole. Motivated by recent X-ray 
cluster studies, we investigate the possibility that local structures at the 100 /i~^Mpc 
scale could be responsible for such correlations. These structures give rise to a local 
Rees-Sciama contribution to the microwave sky that may amount to 1ST jT ^ 10~^ 
at the largest angular scales. We model local structures by a spherical overdensity 
(Lemaitre-Tolman-Bondi model) and assume that the Local Group is falling toward 
the centre. We superimpose the local Rees-Sciama effect on a statistically isotropic, 
gaussian sky. As expected, we find alignments among low multipoles, but a closer look 
reveals that they do not agree with the type of correlations revealed by the data. 
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The microwave sky has presented some surprises at the 
largest angular scales. The Wilkinson Microwave Anisotropy 
Probe (WMAP) confirmed the vanishing of the angular two- 
point correlation function above 60° (Bennett ct al. 2003a), 
a result first obtained by the Cosmic Background Explorer's 
Differ ential Microwave Ra diometer (COBE-DMR) experi- 
ment jHinshaw et al.lfl99d) . In terms of the angular power 
spectrum this implies that the quadrupole and octopole are 
below the theoretical expectation. 

The analysis of full-sky maps that have 
been cleaned from foreground jBennett et alJ 
l2003b^.lFRgmark■ de Oliveira-Costa fc Hamiltonll2003l) . 
lEiriksei^ta7 200^nias revealed further surprises. It was 
pointed out by de Oliv eira-Costa et al. (2004) that the oc- 
topole seems to be planar (all minima and maxima are close 
to a great circle on the sky) and the planes o f the octopole 
and th e quadrupole are closely aligned. lEriksen et alJ 
showed that the northern galactic hemisphere 
lacks power compared with th e southern hemisphere. By 
means of multipole vecto rs llCopi. Huterer fc StarkmanI 
|2004) . ISchwarz et al.l i2004) showed that the quadrupole 
and octopole are correlated with each other and with 
the orientation and motion of the Solar system. The four 
cross products of the quadrupole and octopole vectors are 
unexpectedly close to the ecliptic [> 98% confidence level 
(C.L.)] as well as to the equinox (EQX) and microwave 
dipole (both > 99.7% C.L.) iCopi et aL.2005i) . Based on 
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the additional alignment of a nodal line with the ecliptic 
and an ecliptic north -south asymmetry of the quadrupole 
plus octopole map, ICopi et alJ (|200^ argued that the 
correlation with the ecliptic is unlikely at the > 99.9% 
C.L. In contrast to an unknown Solar system effect, it also 
seems possible that the large-scale anomalies are due to 
a physical correlation with the dipole, in which case the 
correlation with the ecliptic and the EQX would be due to 
the accidental closeness of the dipole and the EQX. 

In this letter we explore the possibility that the ef- 
fect of local non-linear structures on the cosmic microwave 
background (CMB), the local Rees-Sciama (RS) effect 
iRees fc Sciam a' 1968), could induce a correlation between 
the dipole and higher multipoles. In the non-linear regime 
of large-scale structure formation the gravitational potential 
changes with time, and photons climb out of a potential well 
slightly different from the one they fell into. As the CMB 
dipole is considered to be due to our motion with respect to 
the CMB rest frame, and this motion is due to the gravita- 
tional pull of local structures, these structures are a natural 
candidate for contributions to the higher multipoles corre- 
lated with the dipole. Earlier work on the effects of large 
nearby structures on the CMB includes the movi ng cluster 
of galaxies effect JValell2005l : fCoorav fc Setdl2005l) . which is 
different from the effect we are considering here. Second or- 
der corrections to the integrated Sachs- Wolfe effect, which 
catch some aspects of the local RS effect have been consid- 
ered by Tomita (2005a, 2005b). 

The RS effect of distant clusters was estimated to be 
at most 10"'' in a matter-dominated Universe in ISeliakI 



'l99Ci), one order of magnitude below the intrinsic CMB 
anisotropy. The effect of local large structures has been 
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estimat ed to be at most 10~^ usinR t he Swiss Cheese 
model l|Martmez-Gonzalez fc Sand |l990l) and, more reli- 
ably, the Lemaitre-Tolman-Bondi (LTB) model, which is 
the general sph erically symmetric dust solution of the Ein- 
stein equation lfPanek)ll99l ISaez. Arnau fc Fullanal Il993t 
t^lana. Sacz & Arnau 1993i . For an overview. see lKrasMskl 
ltl997D . 

At the time these studies were made, it was gener- 
ally thought that the dipole is mostly due to the infall 
of the Local Gr oup (LG) of galax ies towards the Great 
Attractor (GA) l|Lvnden-BeIl et alJ il988: Drcssleil Il983) . 
a density concentration located 40-60 /i~^Mpc from us, 
with a subdominant component due to the nearby Virgo 
cluster, about 10 /i~^Mpc away. Recent observations of 
X-ray clusters suggest instead that there is a major con- 
tribution to the dipole from the Shapley Supercluster 
(SSC) and other density concentrations at a distance of 
around 130-180 h'^ Mdc JKocevski. MuUis fc Ebeline | 200~ 



Hudson et al]|2004l: iLucev. Radburn-Smith fc Hudsonl200 



iLucev 



.nl2004 



Kocevski fc Ebelinel200.'Tf) . The SSC alone has a density con- 
trast of « 5 over a 30 ft-^Mpc region |Proust ct al. 2005), 
which is 2-3 times the size of the core (of similar density) in 
the GA models. 

An SSC-like object could induce anisotropics at 
the 10~^ level, consi stent with an early estimate in 
iMartfnez-Gonzalez fc Sanz (.l^^( Jp. This can be understood 
by the approximate scaling l|Panekl Il992h 



AT/T ' 



(1) 



where S is the density contrast of the structure, d is its size 
and t is the time at which the CMB photons crossed it. 
For a large angular scale of the source (local and nearby 
structures), this induces contributions to low-^ multipoles, 
especially the dipole, quadrupole and octopole. This could 
include a non-Doppler contribution to the dipole. This would 
imply a change of a few percent in the inferred dipole velo- 
city, which migh t also explain some of the CMB anomalies 
eeman ot al. 200^). The SSC is a non-linear structure, 
and the amplitude of the induced anisotropies cannot be re- 
liably calculated in linear perturbation theory. According to 
a comparison of linear and exac t calculations for GA -like 
objects with the LTB model in iFullana et alJ (^^, lin- 
ear theory is reliable at distances comparable to the Hubble 
scale, but fails for structures within 1000 h~^Mpc or so. 

The advantage of the spherical symmetry of the LTB 
model is that it allows exact calculations for non-linear ob- 
jects; the drawback is that the observed non-linear objects 
such as the GA and SSC do not appear to be spherically 
symmetric. However, one would expect the result to be cor- 
rect within an order of magnitude , and the core of th e SSC 
does seem to be roughly spherical jProust et alj2005h . Also, 
if the preferred direction indicated by the low-^ anomalies 
is due to local structures, this implies that there indeed is a 
degree of symmetry in the local mass distribution. 

In addition, there is a second motivation for study- 
ing a spherically symmetric inhomogeneous model, namely 
dark energy. If interpreted in the framework of isotropic 
and homogeneous cosmology, observations of SNIa imply 
that the expansion of the Universe is accelerating. How- 
ever, in an inhomogeneous spacetime the observations are 
not necessarily inconsistent with deceleration. In particu- 
lar, in the LTB model the parameter go defined with the 



Table 1. Directions of local motion with respect to the CMB 
rest frame. Estimated error for the corrected LG direction of 
IPlionis KolokotroniJ <199Sft rPK1 is 14°, and is 5% for the ve- 
locities. 



Direction 


Galactic coordinates 


V [km/s] 


WMAP dipole velocity 


I = 263°85 ±0°10 


(368±2) 


[Bennett et al. ('20n3a11 


b = 48°25 ± 0°04 




LG velocity 


I = 276° ± 3° 


(627±22) 


[Koeut et al. ("1993)1 


6 = 30° ± 3° 




Virgo infall of LG 


I = 283°92 


170 


[PK (1998)] 


b = 74°51 




Virgo corrected LG vel. 


I = 276° 


510 


[PK (1998)] 


b = 16° 




Shapley concentration 


I = 306°44 




fEinasto et al. C 1997')1 


b = 29°71 





luminosi ty distance is no longer a direct measure of accel- 
eration iHumphrevs. Maartens fc Matraverj Il997^ . It has 
been suggested by several groups that a spherically symmet- 
ric inhomogcnei ty could be u s ed to explain the SNIa data 
( Ce lerierl Il99l iTomital 1200 ll : lAlnes. Amarzguioui fc GronI 
l2005i) . though it is not clear whether such a model could 
be consistent with what is k nown about structures in the 
local Universe llBoleikdl2005li or the observation of baryon 
oscillations in the matter power spectrum. Here we concern 
ourselves only with the CMB. 

The picture of the local Universe that we adopt is 
a spherically symmetric density distribution, with the LG 
falling towards the core of the overdensity at the centre. 
The line between our location and the centre defines a pre- 
ferred direction z, which in the present case corresponds to 
the direction of the dipole after subtracting our motion with 
respect to the LG and the LG's infall towards the nearby 
Virgo cluster (assuming the primordial component of the 
dipole to be negligible). The directions on the sky that are 
important for our analysis are given in Tabled This setup 
exhibits rotational symmetry w.r.t. the axis z (neglecting 
transverse components of our motion). Consequently, only 
zonal harmonics (m = in the z-frame) are generated. 

The density field has two effects on the CMB seen by 
an off-centre observer. First, photons coming from different 
directions travel different routes through the local overden- 
sity, and this creates anisotropy (even with a perfectly ho- 
mogeneous distribution of photons). In a stationary setup 
(i.e. for virialised structures) this effect vanishes and there 
is no imprint on the CMB. Second, the environment will 
affect the evolution of the intrinsic anisotropies (as the ho- 
mogeneous background space does, by changing the angular 
diameter distance). The correct calculation that would ac- 
count for both of these effects would be to study the evo- 
lution of the CMB anisotropies as they travel across the 
density field using perturbation theory on the LTB back- 
ground. We will present the calculation of the amplitude 
of the anisotropies induced by a local structure described 
with the LTB model elsewhere. As in earlier treatments, 
we neglect the second effect and simply add the anisotropy 
generated by the LTB model on top of the intrinsic contri- 
bution. It is possible that this treatment misses some effects 



© 2006 RAS, MNRAS 000,E]ini 



Microwave sky and the local ReesSciama effect 3 




of processing the anisotropies already present. In particular, 
simply linearly adding a new source of anisotropy will in 
general add multipole power, not reduce it, while a proper 
analysis of the processing of the intrinsic anisotropies could 
lead to a multiplicative modifi cation of the amplitu des of 
the low multipoles, discussed in lGordon et alJ (|200^. 

Any other effect with axial symmetry would also in- 
duce anisotropy only for the m — components. It 
has been suggested that spherically symmetric inhomo- 
geneities of the order of horizon siz e or larger would con- 
tribu t e to the low CMB multip oles (|GrishuclLj£.Zerd0vic3 
| l978l: iRaine fc Thoma^ Il98ll : iPacz vhski fc PiranI llQQoF 
iLanglois fc PiranllQQ^) : it was claimed in MoffatI J2005l) that 
this could explain the observed preferred axis. Leaving aside 
the issue that assuming spherical symmetry for the entire 
Universe seems questionable, the observational signature on 
the low multipoles is identical to that from the LTB model 
used to describe local structures, possibly apart from the 
amplitude. 

We study how the CMB is affected by the anisotropy in- 
duced by the additional axisymmetric contribution by means 
of Monte Carlo (MC) maps. We keep the amplitudes as free 
parameters to cover both local and horizon-sized structures 
at the same time, and look at how the observational signa- 
ture compares with what is actually seen on the sky. 

The angular power spectrum is given by Ce ~ 
X^m |fl^m|^/(2^ + l)i where agrn are the coefficients in the 
harmonic decomposition of the temperature map. As pre- 
dicted by the simplest inflationary models, we assume that 
the atm are statistically independent (isotropy), gaussian 
and have zero mean. The aim are then fully characterised 
by angular power. We use the values C2 = 1278.9 /xK^, 
C3 = 590.9 fiK^ from the best-fitting temperature spec- 
trum of the A cold dark matter (ACDM) model with a 
power-law primordial perturbation spectrum to the com- 
bined WMAP, Cosmic Background Imager (CBI) and Ar- 
cminute Cosmology Bolometer Array Receiver (ACBAR) 
dataset.§ For the statistical analysis we generate 10^ MC 
realizations of the quadrupole and the octopole. As the con- 
tribution of the structure described by the LTB model, we 

§ WMAP data products at |http: / /lambda.gsfc.nasa.gov/] 



add to the quadrupole and the octopole a component, de- 
noted by a^o'^') which is a pure m = mode with respect to 
a given physical direction z. 

First we address the amplitude of the quadrupole 
and the octopole. The values of C2 and C3 d eter- 
mined from the WMAP cut-sky jHinshaw et alJlioO^) . the 
TOH map llTeemark et alJl2003D . the Lagrange ILC map 
jEriksen et alj|2004bl) and the Internal Linear Combination 
(ILCl map jBennett et"a r"2003b1 are listed in Tabled The 
extracted quadrupoles have b een Doppler-corrected as de- 
scribed in lSchwarz et al.N2004l) . except for the cut-sky value. 
The values of C2 and C3 from the full-sky maps are signifi- 
cantly larger than the cut-sky values. 

Figure shows how the C2 and C3 histograms compare 
with the data as alQ"*' is increased. For a^J"'' = 40 fiK, 
the number of MC hits that are consistent with the WMAP 
cut-sky data is smaller by a factor of ~ 2 for both C2 and 
C3 compared with the pure CMB sky. For a^o^'^^ = 70 fiK, 
the number of consistent MC hits for C2(C3) is reduced 
by a factor of ~ 5(15) compared with the pure CMB sky. 
Note that adding any power to the theoretically expected 
quadrupole is excluded at the > 99%C.L. level from the cut- 
sky analysis, but for the octopole the same exclusion level 
is not reached until afo"'' = 80 ^K. 

Next we ask what kind of directional patterns the con- 
tribution a^o'^' induces on the s ky. In the multipole vector 
representation dCopi et al any real multipole Tg on a 

sphere (with radial unit vector e) can be expressed with i 
unit vectors v*-*'*' and one scalar A^^^ as 

I I 

7n—~£ i — 1 

The signs of the multipole vectors can be absorbed into the 
scalar quantity A''^\ and are thus unphysical. Note that in 
J^l the r.h.s. contains contributions with angular momen- 
tum (. — 2, £ — 4, ... The uniqueness of the multipole vectors 
is ensured by removing these terms by taking the appro - 
priate traceless symmetric combination jCopi et alJ 1200411 . 
Note that the mult ipole vectors a r e inde pendent of the an- 
gular power. As in 'Sch warz et al.) (|20^), we introduce the 
£{£-l)/2 oriented areas w'*'"'^^ = v'*''' x v^'^'-''. Alignment 
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of the normals n'*^*'-'' = ztw'*'*'-'' with a given direction x 
is tested with 



Sr. 



(3) 



This statistic is a sum over all dot products for a given x, 
so it does not imply any ordering between the terms and 
is a unique and compact quantity. For computi ng the mul- 
tipole vectors we use the method introduced bv lCopi et alJ 
l|2004) .1T 

We look for alignment with three different directions x: 
the north ecliptic pole (NEP), the EQX and the north galac- 
tic pole (NGP). The first two are preferred directions in the 
Solar system and the last defines the plane of the dominant 
foreground. The observed S- values from the different CMB 
maps are given in Table |5| The results of the correlation 
analysis are shown in Fig. |5| By chance the CMB dipole 
and EQX lie very close to each other, so an alignment test 
with the dipole would give results very similar to the one 
with the EQX. 

In the first row of Fig.|21the preferred axis z is chos en to 
be the measured WMAP dipole jBennett et alJl2003a^ . For 
all three tests the anomaly gets clearly worse, i.e. the axial 
mechanism drives the histograms away from the data. Next, 
instead of using t he motion of the L G with respect to the 
CMB rest frame llKoeut et al.lll993l) as the test direction, 
we take the velocity of the LG when corre cted for Virgocen- 
tric motion lonis fc Kolokotronij|l998ft . since this differs 
more from the WMAP dipole. The results are shown in the 
second row of Fig. |21 The situation for the alignment with 
the EQX is again worse, but there is not much effect on the 
ecliptic alignment. For the alignment with the galactic plane, 
the axial contribution makes an apparent Galactic correla- 
tion more probable, i.e. there is a certain probability of over- 
estimating the galactic foreground. For both test directions 
(rows one and two), the alignment with the EQX gets worse. 
For example, in the direction of the Virgo-corrected LG mo- 
tion an exclusion of ~ 99.9% C.L. for can 
be given for all three maps. Note that adding any multipole 
power in this test can already be excluded at the ^ 99.4% 
C.L. 

As a complementary test we show the alignment likeli- 
hood with regard to an orthogonal test direction, the NEP, 
in row three of Fig. |5| An ecliptic extra contribution in the 
CMB would indeed induce an alignment of normal vectors 
similar to the observed one. In particular, for aJo'^' ~ 50 mK, 
the probability of finding an alignment with the NEP itself 
becomes roughly 5%, and the probability for the EQX align- 
ment rises to 1%. 

To summarise, the results of recent X-ray cluster studies 
indicate the existence of large-scale fluctuations S ~ 0(1) at 
distances of ~ 100 /i~^Mpc. The local RS effect on the pri- 
mordial photons from these structures is estimated to be of 
order ~ 10~^ at large angular scales. This raises the question 
of whether the local RS effect can account for the observed 
anomalies in the low multipoles of the CMB. 

In this letter we have assumed spherical symmetry of 
the local superstructure, with an object like the Shapley Su- 
percluster at the centre. Under this assumption we should 
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Table 2. Tests, defined in j3J and explained in the text, appUed 
to the TOH, Lagrange ILC and ILC maps. All quadrupoles except 
the cut-sky value have been Doppler-corrected. 

Cut sky TOH map LILC map ILC map 



SnNEP 
•S'nEQX 
S'nNGP 



129 AtK^ 
320 ^iK^ 



203 AtK^ 
454 ^iK^ 
0.194 
0.886 
0.803 



352 AtK^ 
571 AtK2 
0.193 
0.866 
0.803 



196 iiK'^ 
552 fiK^ 
0.210 
0.870 
0.810 



observe an axisymmetric effect on the microwave sky. The 
preferred axis has been taken to point in the direction of the 
local velocity flow (not shown in Fig.|5J, the CMB dipole and 
the Virgo-corrected LG flow vector. We have added this ax- 
isymmetric contribution to a statistically isotropic gaussian 
random map and compared it by means of the S-statistic 
with WMAP measurements. The result is that this mecha- 
nism can be excluded at > 99%C.L. Our analysis also applies 
to any other effect which gives an axisymmetric addition to 
the statistically isotropic and gaussian random sky. How- 
ever, we find that a Solar system effect would be co nsistent 
with the data, in agreement with lCopi et al.l (120051) . 

The present work also indicates that the local Rees- 
Sciama effect might be important for the interpretation of 
WMAP data and future PLANCK data on the largest an- 
gular scales. Our work suggests that a more detailed study 
of the Rees-Sciama effect could lead to useful constraints on 
the local large-scale structure of the Universe. 

Supplementary figures, illustrating the RS- 
induced CMB correlations, are available at 
http:/ / www.physik.uni-bielefeld.de / cosmology/rs.html, 
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Figure 2. Alignment statistic JJJ foi' quadrupole and octopole normals. For the three columns we pick the test direction x to be NEP, 
EQX and NGP respectively. In rows we consider three different choices of the preferred axis z, namely the WMAP dipole, the direction of 
motion of the LG, corrected for its Virgo infall, and the direction of the NEP. Shown are the likelihoods of the 5-statistic for statistically 
isotropic gaussian skies (thick solid lines) as well as different magnitudes of axial contamination of the CMB. Vertical lines represent 
the measured 5-values from the TOH (solid line), LILC (dotted line) and ILC (dashed line) maps (see Table 1^. Introducing a preferred 
axis induces correlations. For the directions of local motion (first and second rows) these correlations make the discrepancy between the 
measured S-values and model even bigger. At the same time, a Solar system effect is more consistent with data. 
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